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Abstract 

We analyse how a probing particle modifies infrared electromagnetic near fields. The particle, 
assimilated to both electric and magnetic dipoles, represents the tip of an apertureless scanning 
optical near- field microscope (SNOM). We show that the interaction can be accounted for by 
ascribing to the particle effective dipole polarizabilities that add the effect of retardation to the 
one of the image dipole. Apart from these polarizabilities, the SNOM signal expression depends 
only on the fields without tip perturbation, shown to be closely related to the electromagnetic 
density of states (EM-LDOS) and essentially linked to the sample's optical properties, so that 
measuring local spectra of heated samples is equivalent to performing a local surface spectroscopy. 
We also analyse the case where the probing particle is hotter. We evaluate in this case the impact 
of the effective polarizabilities on the tip-sample near-field radiative heat transfer. We also show 
that such an heated probe above a surface also performs a surface spectroscopy. The calculations 
agree well with available experimental data. 

PACS numbers: 07.79.fc,44.40.+a,71.36+c 
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I. INTRODUCTION 



Since the seminal works of Rytov^, it is known that thermal radiation has a different 
behaviour when the involved caracteristic lengths are large or small compared to the ther- 
mal wavelengthPSl. For example, the heat flux transferred between bodies separated by a 
subwavelength distance can exceed by far the one between perfect blackbodies^. Energy 
density^ and coherence properties^ are also strongly affected in the near-field, especially close 
to materials exhibiting resonances such as polaritons. Knowing precisely how the electro- 
magnetic field behaves close to a surface is therefore an important issue in order to address 
all potential applications involving near-field heat transfer. 

From an experimental point of view, near-field radiation coherence properties have been 
utilized to produce directional and monochromatic thermal sources^^. More recently, near- 
field radiative heat transfer enhancement have been detected by means of probe microscopy 
techniques'^. Near-field thermal flux imaging has been operated with a scanning thermal 
microscop e 1 16 * 17 *. A scanning near-field optical microscopy (SNOM) without external enlight- 
ment, termed thermal radiation scanning tunneling microscopy (TRSTM), has also been 
shown to be useful to image surfaces^ and to be very sensitive 1 19 * 20 *. Very recently, local 
spectra has also been performed^. All these experimental techniques use a small probe 
brought to the vicinity of the studied material surface. When the probe is approached 
close to the surface, mutual interaction between the tip and the surface modifies the local 
electromagnetic fielcP^^, a phenomenon described by probe optical properties such as its 
polarizability if it is modelled as a dipole. This is the origin of the field scattering to the 
far field, which is beneficial for TRSTM, but also complicates the data analysis for the 
previoulsy-mentioned near-field techniques although one can relate the local electromag- 
netic field to the surface material optical properties. A fundamental issue is to know how 
the TRSTM signal scattered by a tip and detected in the far-field can be related to the 
electromagnetic local density of states (LDOS)^^3. Is there a way that a SNOM detecting 
thermal radiation might be the electromagnetic equivalent of the scanning tunneling micro- 
scope detecting electronic LDOS^? Moreover, when a probe is heated, one can ask itself 
how fast it cools down and what the signal detected in far-field is when it is approached 
close to a surface. 

If some of these questions have already been addressed in the past! 26 1 28 1 , our goal is here 
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to clarify remaining interrogations and to complete the body of work. Following previous 
similar works^^, we will first see how the particle polarizability can be replaced by an 
effective polarizability taking into account the multiple reflections between the probe and 
the surface. This model is more general than the image-dipole one^ whose range of validity 
is very restricted in the infrared. We will then use the theory to calculate the signal detected 
in the far-field due to scattering of the near field by a probe. We will simulate the SNOM 
signal obtained by scanning a surface excited either by a plasmon or more generally by 
thermal excitation (near-field thermal emission). We will consider in the following two parts 
that the probe tip can be assimilated to both electric and magnetic dipoles. We will first 
give the signal detected in far field due to a heated probe when accounting for the surface 
reflections. We will show then how cooling is increased when the particle is approached close 
to the surface. 



II. EXPRESSIONS OF THE EFFECTIVE POLARIZABILITIES 

We propose here to calculate the effective polarizability of a dipolar particle when it is 
placed in an environment which is different from vacuum. Indeed, when a particle is added 
in a system, the electromagnetic field present in the system illuminates the particle, what 
gives birth to an induced dipole moment centered at its position. This dipole radiates also 
a field everywhere, including at the particle position. Interactions between the particle and 
the system modify the total electromagnetic field, which becomes then different from the 
one in absence of a perturbating particle. Our aim is to show that we can consider the 
electromagnetic field to be unchanged at the particle's center if we ascribe to the particle 
an effective polarizability that accounts for the perturbation (Fig. [I]). 

Let us call E° the electromagnetic field in the system in the absence of particle. When a 
particle is placed in the system, an electric dipole p and a magnetic dipole m will appear 
at the position of the particle r p . These dipoles radiate a field : the total local field E' oc is 
the sum of the field in the absence of particle with the field radiated by the dipoles. In the 
following, we write the dipole radiated field with the Green's tensor help. 

E' oc (r) = E°(r) + G^(r, r p ) • p + G^(r, r p ) ■ m (1) 
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Particle 




FIG. 1. System scheme 

that can also be written 

E loc (r) = E°(r) + G^(r, r p ) ■ a E E loc (r p ) + (f" (r, r p ) ■ a H H loc (r p ) (2) 

where a E and a H are the particle electric and magnetic dipoles when it is alone in vacuum. 
Analogous expressions exist for the magnetic field 

4* HE o HH 

H (r) = H (r) + G (r, r p ) • p + G (r,r t )-m (3) 

that can also be written 

H ioc (r) = H°(r) + G^(r, r p ) • a E E loc {v p ) + G^(r, r„) • a H H loc (r p ) (4) 

The Green tensors used here are of four type : G (r, r p ) is the Green tensor that gives 
the electric field at position r when an electric dipole is placed at r p . In the same way, 
G (r, r p ) is the Green tensor that gives the electric field at position r when an magnetic 
dipole is placed at r p . G and G respectively gives the magnetic field of an electric and 
a magnetic dipole. These Green tensors, which can all be expressed in terms of G (r, r p /22, 
are the sum of a direct contribution (that is the Green tensor in vacuum) and of a perturbed 
contribution. When one considers the electromagnetic field at the particle position r p , only 
the perturbed Green tensor contributes. Then, the local field is obtained solving the system 




E te (r p ) = A _1 E°(r p ) + B H°(r p ) (5) 
H ioc (r 4 ) = C _1 E°(r p ) + D^H (r p ) (6) 
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where 



A=[l-a E G (r p , r p )] - a H G (r p , r p )[l - a H G (r^r^^/G (r p ,r p ) (7) 

/ ) j-^HE O <r>HH ir^EH 4-i ; > -EE 1 

B = -a £ G (r p ,r p ) + [I -a H G (r p , r p )][c^G (r p , ip)]" 1 ^ - « S G (r p , r p )] (8) 

c 



-c^G 



r p ,r p ) + [l -a £ G (r p ,r p )][a £ G (r p , r,)]"^! - a H G (r„,r p )] (9) 



D = [l-a^G (r p ,r p )] -« E G (r p ,r p )[l 



„^EE rr^EH 

a E G {r p ,r p )]- l a H G (r p ,r p ) 



(10) 



The preceding equations are general for the calculation of a local field at dipole position. 
Apart from the fact that magnetic dipole is taken into account, the self-consistent reasoning 
used to obtain local fields is very similar to previous work using the so-called coupled dipole 
theory such as ^ 24 * 32 l Local fields only depend on the system Green's tensors. Particle 
dipoles can be related to the field at the tip position in the absence of particle 

. — , . — \ / 

r, i 

(11) 




' a E A a H B 



a E C a H B 




so that the four elements in the matrix can be seen as effective polarizability expressions. 

We now highlight the case of a simple system made of a planar interface separating a 
material and vacuum. In this situation, Green tensors are well knowrpE2] anc j can De written 
knowing the materials optical properties. We recall the direct Green's tensor describing 
propagation in vacuum 



► AB 



G (r,r,w) 



iklC AB 



d 2 K 



e A e B 



57T Z J 7 

and the one associated to a plane wave reflection at an interface 

d 2 K 



e P- e P-yK.(R-R') e -n(z-^') 



++AB 

G R (r,r,w) 



iklC AB 



e s A + r s e s B - + e^' r' J e 



P+ r P e P- \ „iK.(R-R') J-y(z+z') 



B 



(12) 



(13) 



»7r' J 7 

Here K and 7 are respectively the parallel and the perpendicular component of the wavevec- 
tor. r s and r p are the Fresnel reflection coefficients for polarizations s and p which depend 
only on the material optical properties and K. Introducing vectors s and p~ vectors used 
by SiptP : 



s = K x z/K 



V 



k K 
1 



K 2 z + jK x ± + jK y y 
K 2 z - jK x ± - 7-fCy 



(14) 
(15) 

(16) 
(17) 



one writes the vectors in Green dyadics as 



P P ~ 



E 



-p 
-P 



(18) 
(19) 
(20) 



Moreover, Cee = £o 5 Ceh = Hoc, Che = eoc, Chh = 1- When one considers the elec- 
tromagnetic field at the particle position r p , only the reflected Green's dyadic contributes. 

^EH ^HE 

In this plane-parallel geometry, G R (r p , r p ) = G R (r p , r p ) = 0. Expressions of the electro- 
magnetic field at the tip position can be obtained from simplification of (|5| and ^ taken 



in r = r„. Then 



H /oc ( 



<-> F ^EE 

I - a G R (r p ,r p ; 

O HH 



I - a H G 



R 



v r p ; r pj 



E°(r, 



H°(r p ) 



(21) 
(22) 



leading to the following dipole expressions : 



P(r P 



m r r 



a 



,^EE 



a 



H 



I - a G R (r p ,r p ) 
I-a G R (r p , r p ) 



E u (r r 



H°(r r . 



(23) 
(24) 



We note that we can recast these formulae in order to directly relate the dipole moments 
and the fields in absence of dipole by introducing the following effective polarisabilities 



a 



a 



I -a E G EE (r t ,r t ) 



(25) 



and 



a = a 



I - a H G HH {r t ,r 



(26) 



These effective polarizabilities are anisotropic and account for the fields multiple reflections 
between the particle and the interface. 



III. PARAMETRIC STUDY OF THE EFFECTIVE POLARIZABILITIES 



In this section, we analyze how the effective polarizabilties vary as a function of the ma- 
terial optical properties, particle sizes and particle-surface distance. We consider that the 
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particles can be assimilated to spherical dipoles, with well-established "vacuum" polarizabil- 
ities. Note that these ones depend, independently from the influence of the interface studied 
here, on the ratio of the particle radius R p to the wavelength A (Mie parameters = 2tcR p / A) 
and also on y = yfex, where e is the particle dielectric constant. The following expressions 
given by Chapuis 34 apply if the wavelength is much larger than R p : 



a E (uj) = e^ixRl 



2 [sin(y) — ycos{y)] — x 2 


— sinfw) i cos(v) i • / \ 

[ y 2 ' + y ' + sm(y)\ 


[sin(?/) — ycos(?/)] + x 2 


— sin(-y) i cos(i/) i • / \ 
y 2 + y + sin ( J/) J 



and 



a H (u) 



-2ttRI 



1 - 



x 
10 



+ 



3 3 

+ -cot(y) 

y y 



l 



X 



(27) 



(28) 



Our aim is to compare the theoretical results to apertureless Scanning Near-field Mi- 
croscopy (SNOM) results. We consider particles made of tugnsten, as this material is often 
used for SNOM tips. We study the effective polarisabilities for two spherical tip sizes (lOOnm 
and 1 /zm radii) and above three materials : SiC and Si02, being both dielectrics, and gold, 
a metal. 



A. Different contributions to the effective polarizabilities 

Let us consider the case of a lOOnm-radius tip of tungsten situated 100 nm above an 
interface separating the vacuum and the sample. This corresponds to a sphere in contact 
with the interface. In principle, the dipole model is not sufficient anymore to describe 
this situation since the evanescent field is not constant inside the particle and therefore 
multipolar contributions should be considered. However, due to the complexity of this task, 
the use of dipoles to model the tip has been the rule for several SNOM theories and will be 
considered as a limiting case here. It shows most of the modifications that occur when the 
tip is approached close to a surface. 

In the case of a dipolar particle close to a plane interface, four different contributions to 
the polarizabilities are identified: The parallel electric (a^ x ), perpendicular electric (af 2 ), 
parallel magnetic (a^ x ) and perpendicular magnetic (af z ) ones. They are plotted in Fig. [2J 

We observe that both perpendicular and parallel effective electrical polarisabilities are 
very different from the vacuum polarisability of the single particle. In particular, resonances 
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FIG. 2. Parallel electric (a), perpendicular electric (b), parallel magnetic (c) and perpendicular 
magnetic (d) polarizability for a lOOnm-radius spherical tungsten tip situated at z=100 nm above 
an plane suface of SiC, Si02 or Gold. 

are observed respectively close to the SiC and Si02 plasmon polariton resonance frequencies. 
This is not surprising since the single-interface Green tensors enter in the expresions of the 
effective polarizabilities. Since the dyadics depend strongly on the reflection coefficients, 
the former take large values close to the resonance frequencies, especially at small distances 
from the interface. Thus the effective polarizabilities are greatly affected at small distances 
and close to resonant frequencies. In contrast, no peak is seen in the spectra of gold effective 
polarizabilities since gold does not exhibit resonances in the studied frequency range. Thus, 
the effective polarizabilities are different from the vacuum ones as the presence of evanescent 
modes close to the interface affect the particle. The observation of the effective magnetic 
polarizaties spectra show very different variations from the electric one. The main reason is 
that vacuum magnetic polarizabilities have different asymptotic variations at small radii R p . 
Indeed, the electrical polarizabilities behave in the small-particle regime as Bl (Clausius- 
Mossotti limit) whereas the magnetic ones behave as Rp k$ (e — In the present case, 
as the particle radius is much smaller than the wavelength, the magnetic polarizability is 
smaller, which implies that the correction is also smaller. We note that the correction is 
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FIG. 3. Parallel electric (a), perpendicular electric (b), parallel magnetic (c) and perpendicular 
magnetic (d) polarizability for a lOOnm-radius spherical tungsten tip situated at z=300 nm above 
an plane suface of SiC, Si02 or Gold. 



larger for gold than for SiC and SiC^. This is a consequence of the fact that evanescent 
magnetic modes are more present above metals than above dielectrics, as shown by the large 
value of the metal dielectric constant in the infrared. 

The same particles polarizabilities have been studied at a larger distance to the interface 
i.e for 300 nm. The variations are shown in Fig. [3] The electrical polarizabilities spectrum 
behaviour with variations around resonant frequencies is still present at larger distance 
although the corrections to the vacuum polarizability are much smaller. Thus, evanescent 
modes decay exponentially when the distance to the interface increases. This fact also 
explains why corrections are almost inexistant at z=300 nm for magnetic polarizability. 

Let us now study the case of a larger particle (1 /mi) situated at 1 /im form the interface. 
The effective polarizabilities are represented in Fig. |4j We note important corrections 
for both effective electric and magnetic polarizabilities. These corrections are once again 
much more important around surface resonance frequencies. The impact on the magnetic 
polarizabilities is now much more striking than for small particles due to the fact that here 
k R p is larger. We will see in the next section that retardation effects begin to enter into 

10 




FIG. 4. Parallel electric (a), perpendicular electric (b), parallel magnetic (c) and perpendicular 
magnetic (d) polarizability for a 1/xm-radius spherical tungsten tip situated at z=l fim above an 
plane suface of SiC, Si02 or Gold. 



play at such distances so that corrections at 1/iin are different from the one at 100 nm as 
can be noticed when one analyzing carefully the effective polarizabilities correction curves. 



B. Effective polarizability asymptotic expression 

We now study the effetive polarizabilities in order to obtain asymptotic expressions in 
some limited physical situations. We first start with Green tensors that appear as an integral 
over the parallel wavevector: 



<-> EE 

G R (r t ,r t ) 



roo KdK ( s 
J 2-y v 



47T 



,p7 \ p 2i^z t 
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oo KdK / s 
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oo K z dK 



IK 



(29) 
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and 



G R (r t ,r t ) 



( roo KdK/ p _ 
JO 27 V 



Aire 2 



\ 








roo KdK , 








•00 K 3 dK s,2i"/z t 
7 fcg r e 



(30) 



We suppose now that the separation distance is much smaller than the wavelength, so 
that the distance-dependent terms can be replaced by their asymptotic expressions. In that 
case, r p = (e — l)/(e+ 1) and r s = (e — V)k\jAK 2 . Integration over the parallel wavevector 
is easy since r y 2 = k 2 — K 2 ~ — K 2 so that 7 = iK. We recall that this is valid in a regime 
sometimes called "extreme near-field"^, where the integrals are mostly due to their large K 
contribution, which is not always the case®' even in the near field. In this case, asymptotic 
expressions of the Green dyadics read 



^EE 



32ttzI 
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So that the effective polarizabilities take the form 
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(34) 



These expressions are the asymptotic approximations of the Green tensors. For the electric 
term, the obtained expression is also known as the electrostatic limit, which means that 
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retardation is not taken into account. Note that the magnetic terms vanishes if retardation 
is discarded due to a H . 

The resulting expressions for the electric effective polarizabilities are very close to previous 
work of Knoll and KeilmamP^. In Fig. [HJ the example of a tungsten spherical particle above 
SiC>2 shows that this approximation is very good for small particles at short distances. For 
larger particle, such as a 1/xm particle located at a 1 fim distance, some deviation exists 
between the exact calculation and asymptotic expressions. This is due to retardation effects 
that enter into action and which are clearly not negligible even at subwavelegth distances 
as shown here. 

Fig. [5] shows also that the effective magnetic polarisabilities and their approximations fit 
together at shortest distances but not anymore at larger ones. However, we note that for 
small particles, the correction at large distances (but still in the near field) is negligible due 
to the smallness of the polarizability. We underline also that the electrostatic approximation 
neglecting any magnetic term can be used only at extremely small distances. 

IV. APPLICATION TO SNOM DETECTION 

The spherical particles precedently studied model tips in SNOM experiments. The goal 
of this kind of experiments is to study the electromagnetic field close to a surface, the tip- 
surface distances being much smaller than the wavelength considered as a consequence (near 
field). The tip scatters the near-field electromagnetic radiation, which is then detected in 
the far-field by a detector located at a position r^. The detection scheme is represented in 
Fig. § 

A. Signal at the detector 

The tip possesses electric and magnetic polarizabilities, so that the incoming electromag- 
netic field creates electric and a magnetic dipole moments at the tip position r t . These 
dipoles radiate an EM field which is detected at the detector position, in the far field. Thus, 
the detected field is just the simple electromagnetic field radiated in far field by an electric 
and a magnetic dipoles taking into account the influence of the interface. The signal at the 
detector is calculated in the so-called far-field approximation, where the distance has to be 
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FIG. 5. Comparison between effective polarizabilities obtained with the exact theory developed 
here and obtained through electrostatic hypothesis. Four top figures : lOOnm-radius W sphere 
situated at z = 100 nm of a surface of Si02. Four bottom figures : 1/xm-radius W sphere situated 
at z=l /jm of a surface of Si02. 
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large compared to square of scatterer size p over the wavelength: d » p 2 /\. The signal at 
the detector can be considered as a plane wave so that it reads 

(S(u>)) = e -°^\E d (u)\ 2 r 2 dn (35) 

Let us introduce, the elementary vector of the tip-detector direction u d = (r d — r t )/\r d — 
r t \ = u d \\ + u d ± where u d ± is perpendicular to the interface. We define u d = u d \\ — u d ±, 
s d = u d \\ x e z , = u d x s d and = — x s d . Let us also introduce four more tensors 
h(ud), g(u d ), h (u d ) and g (u d ). 

h(u d ) = sJd + pjpj- (36) 

g(u d ) = -s d pj + pjs d (37) 

h (u d ) = (s d r s d s d + pjr p d p^ (38) 
++R 

g \ Ud ) = (-s d r s d p d +pjr p d s d )e^ (39) 

where denotes the phase difference between radiation going from the tip to the detector 
and the one going from the tip to the detector after one reflection on the surface. Following 
classical expressions of field radiated by dipoles in the far field, the electromagnetic field 
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radiated at the detector is 
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If we introduce now r (u d ) = h(u d ) + h (u rf ) and r (u d ) = gf 



(40) 



Ud) + g (ud), we find, 



replacing (35) the expression of the power received at the detector in a solid angle dVt around 
the detector direction: 



(S) 



32n 2 c 



(41) 



The final result is a combination of the different electromagnetic field components. These 
expressions are quadratic expressions of the electromagnetic field component and are related 
to the electromagnetic energy or the Poynting vector close to the interface. 



B. Expressions of the tensors h and g. 



We consider now a situation where the detector makes an angle 9 with the vertical axis 
in the z — y plane. Therefore, in spherical coordinates, the detector is at position (r^, <p, 6) 
with ip = it/ 2 (Fig. m). The tip is at position z t above the interface. 



Detector 



X 




(p=rc/2 



FIG. 7. Tip and detector position 
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At first order in zt, the tensors reads 
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(42) 

-sin0(l + r s (0)e 2ifeo2tcose ) \ 


/ 

(43) 



V. APERTURELESS SNOM SIGNAL FOR VARIOUS SITUATIONS OF INTER- 
EST 

A. Signal generated by a surface plasmon excited with a quantum cascade laser 



We now calculate the signal detected by an apertureless SNOM above a material sup- 
porting surface plasmons. Typical materials are metals^ such as gold and all materials for 
which the dielectric constant is smaller than -1. Surface plasmons are eigenmodes of a planar 
surface. They correspond to a pole of the plane interface reflection coefficient. They only 
exist for TM (or p) polarization. If z is the direction perpendicular to the interface, the 
plasmon magnetic field is given by 

H(r) = (44) 





whereas the electric field reads 



E(r) 



70-^0 e i(Ky+~(oz) 
KH i(Ky+-y n z) 



(45) 



If the detector makes an angle with the vertical axis and is situated in the z — y plane, 



the signal at the detector is 
(S(r d )) 



32tt 2 c 2 



(46) 
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In upper Fig. [8j we plot the signal at the detector with the tip-sample distance in the case 
of a 1 /im-radius spherical dipole which simulates a ungsten tip oscillating in the vertical 
direction with an amplitude of 150 nm. These conditions are the ones used in SNOM 
experiments such as Reffi We observe oscillations as the tip is retracted away from the 
surface. This is due to interferences between the signal scaterred by the tip and radiated 
directly to the detector and the signal scattered by the tip that undergoes a reflection before 
reaching the detector. The lower curve in Fig. [8] shows the same kind of signal in the case 
in which the tip is not modulated. The signal behaviours obtained from an oscilating tip 
or an non-oscilating tip are similar. The oscillation period is around A p / asmon /2 = 3.75 /zm 
and the signal decreases with the tip-sample distance. We observe a similar behaviour when 
one probes the plasmonic signal above gold as seen in Fig. [9j Here the plasmon is excited 
by a quantum cascade laser at 7.5 fim. Besides the oscillations, the difference between 
the theoretical and the experimental curves may be due to different geometry of the signal 
collection, which involves a Cassegrain objective in the experiment, and to the fact that 
our model approximates the tip by a 1 m radius tungsten sphere while it is in reality made 
scatterers distributed along a long cone with a curved apex. 



B. Signal due to thermal emission 

1. Thermal fields above a sample 

When the electromagnetic field above the surface is only the one due to thermal emission 
from the sample surface, several simplifications occurs. Cross correlations functions such 
as (E° x (v t )E° y *(v t )), (JS2(r t )^(r t )>, (E° y (r t )E° z *(v t )), (#°(r t )# y °*( r *)}> (H°(r t )H»(r t )), 
(H y (v t )H° z *(r t )), (E° x (r t )H° z *(v t )), (E° z (r t )H^(r t )), (E%T t )H?(r t )) and (E°(r t )H»(r t )) 
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FIG. 8. Theoretical signal above a sample of gold on which plasmons have been excited at a 
wavelength A = 7.5 fj,m. Upper figure : signal detected by a l^m-radius tungsten tip above a 
gold plane interface supporting a plasmon excited at A =7.5 /im. The tip oscillates with a 150 
nm-amplitude. The amplitude of the first harmonic of the signal is plotted, corresponding to the 
signal detected at the tip oscillation frequency. Lower figure : same conditions as above except 
that the tip position is not modulated 

are all equal to due to the fact that thermal currents are decorrelated for different 
directions^. Moreover, due to rotational symmetry around any axis in the z direction, 
(\E° x (rt)\ 2 ) = (\E° y (v t )\ 2 ) and <|tf°(r t )| 2 > = (|^°(r,)| 2 ). 

In these conditions, the signal at the detector in a direction making an angle 6 with the 
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FIG. 9. SNOM signal above a gold sample supporting a plasmon excited by a quantum-cascade 
laser at 7.5 /im. Detecting tip is made of tungsten and oscilation ampitude is 150 nm. 

z axis reads 
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Expressions of the electromagnetic field quadratic quantities due to thermal emission above 
a plane interface are well known in the literature^' : 
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where v = yl — u 2 and Q(uj,T) = l/[exp[fwj / (ki,T)] — 1] is the mean energy of an oscillator 
at angular frequency u at thermal equilibrium. As a consequence, if the substrate and 
the particle dipole materials are known, the signal at the detector coming from an emitted 
thermal electromagnetic field is known. If the signal is divided by the mean energy of an 
oscillator Q(u,T), it only depends on the particle polarizability, the tip-sample distance 
and the surface optical properties. If the polarizability is known, one can map a quantity 
which is specific of the surface and therefore make a surface spectroscopy. Note however that 
expressions of the effective polarizability and of the thermal electromagnetic field are closely 
related. Indeed, both are depending of the system Green tensors taken at the tip position. As 
a consequence, if there is a resonance at a given frequency for the quadratic electromagnetic 
field expressions, it is likely that there will be also one in the effective polarizability spectra, 
at the same frequency. This constitutes a difficulty for the interpretation of the signal 
at the detector. Indeed, quadratic electromagnetic quantities are closely related to the 
electromagnetic LDOS (EM LDOS)Pl Relating the signal detector to the LDOS would be 
very interesting since it would give a way to detect this quantity as electronic tunneling 
microscopy does it for electronic LDOS 27 . Unfortunately, it is not possible to find a simple 
and universal relation between the two quantities. However, in some specific situations, 
when one term is leading, the EM LDOS and the detected signal can be proportional. 



2. Probing a polar material 

a. Signal simulation We now study the spectrum of the signal detected by a probe 
above SiC. As we have seen in the previous section, this signal strongly depends on the 
effective polarizability and therefore on the tip size and the tip distance. To illustrate this, 



we show in Fig. [10] the spectral signal detected by a lOOnm tungsten probe above SiC and 
by a 1/^m-radius probe. It is observed that the emission peak is rather narrow for the small 
tip. This peak is also situated around the phonon polariton resonance frequency that is 943 
cm -1 for SiC. We also note that for a small tungsten tip, the signal above SiC is dominated 
at short distance by electrical terms. We see that both parallel and perpendicular terms 
have a significant contribution. These contributions are slightly shifted with respect to each 
other: Effective parallel and perpendicular resonances are indeed different with a factor of 2 



(Eq. 33 ) in the right terms of the denominator asymptotic expressions, leading to a different 
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FIG. 10. Signal spectrum at the detector. The tip is constituted of a W spherical particle situated 
above SiC. The signal is detected in a direction making an angle of 30° with the vertical (a) Total 
signal and different contributions for a lOOnm radius tip at z =100nm above the surface, (b) Total 
signal and different contributions for a 1 //m radius tip at z =1 /um above the surface, (c) Signal 
spectrum for z = Ru p = 100 nm, z = Ru p = 1 //m and z = Ru p = 1.6 /im. Comparison with 
LDOS (d) Signal spectrum for Ru p = 100 nm and two different distances to the interface z = 100 
nm and z = 1.6 ^m. 



frequencies for the resonant peak 

For a 1 fim tip situated at 1 /im from the interface, we note that the main contribution 
comes from the parallel electric term. Magnetic terms are also more important here than for 
small tip due to the fact that the magnetic polarizability is larger. These magnetic terms 
do not contribute a lot because magnetic energy density is lower than electric energy near 
the polariton resonance for polar materials. However, the mixed term, which involves both 
magnetic and electric terms is here the second contribution to the signal and is not negligible. 



Fig. 10 (c) shows various signal spectra for tips that would be in contact with the surface. 
They are compared with the theoretical spectrum of the LDOS at 100 nm. A lOOnm-radius 
tip gives a signal similar to the electromagnetic LDOS: A well-defined peak appears around 
the polariton frequency. The peak width is very close to the LDOS one but the peak position 
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is nevertheless slightly shifted. When the tip size increases, the peak tends to shift and to 
broaden. Recent SNOM experiments based on the measurement of the near-field thermal 
emission using a tungsten tip seem to confirm this shifiting and broadening suggesting that 
the approximation of the tip by a simple spherical dipole which we make in our model is 
valid to some extent. Note also that this broadening, although less pronounced than for a 
large tip, also occurs when a small tip is retracted from the surface as it can be seen in Fig. 



10 (d). 



b. Relation between the detected signal and the LDOS As shown in Fig. 10, the signal 
calculated with a small tip is very similar to the LDOS spectrum around the polariton 
resonance. At short distance, the signal is indeed mainly dominated by the parallel electric 
contribution. This means that the signal mainly depends on a^ x and on \E®\ 2 . In the case 
of a thermal signal, this last quantity is representative of the electrical energy density and of 
the electromagnetic LDOS. Therefore, a SNOM experiment detecting thermal near field will 
measure the product of one of the effective polarizability terms by a partial contribution to 
the EM LDOS. If the effective polarizability has a flat response i.e. its does not vary with the 
frequency (or only weakly), one can say that the signal at the detector is proportional to the 
EM LDOS. This is not the case in the present situation where polarizability can be increased 
by a factor of 3 to 5 around the resonance: The resulting signal is the product of two peaks 
situated approximately at the same frequency. As the value of the EM LDOS in the peak 
spectral band is abut 100 to 1000 times its value outside this band, its multiplication by 
the effective polarizability will give a peak that is not exactly the EM LDOS but which is 
representative of the LDOS. 

We now consider larger tips of micrometer size radii. The signal at the detector scattered 
by such tips has a broader spectrum and the frequency corresponding to the emission spec- 
trum is shifted to lower frequencies. When the tip size increases, polarizability is shifted 
to lower frequency and broadened. A larger distances, EM LDOS is also broadened and 
shifted but differently than the polarizability. Frequency shift is indeed less important for 
EM LDOS. The resulting signal is now the product of two peaks which are rather similar 
in terms of broadness but with their maxima separated. The resulting signal is a peak situ- 
ated between the polarizability and LDOS peaks. With these conditions, the signal at the 
detector is a peak related to the EM LDOS but which is not stricly-speaking the EM LDOS. 
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3. Probing a metallic sample 



We now consider the example of a tip situated above a metal surface. Fig. 11 shows the 
signal scattered by a tungsten particle above a gold surface heated at 300 K. For a small 
particle, the contribution to the signal is dominated by the magnetic term below 1000 cm -1 . 
The main reason for that is the fact that the magnetic energy is more important than the 
electrical energy in this part of the spectrum where Au is highly reflecting. As seen in Fig. |2j 
the magnetic polarizability has a flat response for small tips above 400 cm -1 . Then, between 
400 and 1000 cm -1 the signal at the detector is proportional to parallel magnetic energy 
which is close to total electromagnetic energy. Thus, the signal detected in this spectral 
range is close to the EM LDOS. On the contrary, above 1000 cm -1 , the signal is dominated 
by the perpendicular electric contribution but the parallel magnetic and the mixed term 
are of the same order of magnitude even if smaller. For such situation, it is not correct to 
state that the signal is proportional to the EM LDOS although it is related to it. For a 
large tip, there is no spectral range where a term clearly dominates the contribution to the 
signal (except at high frequency where the validity of the dipole approximation becomes 
questionable). The detected signal cannot be considered as simply proportional to the EM 
LDOS. 



VI. THERMAL EMISSION OF A NANOPARTICLE 

We now consider a small dipolar particle heated at temperature T. This particle radiates 
an electromagnetic field in all space. This field is a thermal field which is here detected in 
far-field by a detector. The signal at the detector has a contribution coming directly from 
the particle in straight line or after one reflection at the interface. The signal at the detector 
has the following expression: 



(S(u)) 



Q(co,T)^ [S[«fj (|1 + r ^e)e 2ikoztcose \ 2 + cos 2 9\1 - r p (9)e 



P( n\ a 2ikg zt cos 8 1 2 



327T 3 ' ' c 

+ 3 [afj 1 1 + r p (6)e 2ikoZt cos9 \ 2 sin 2 6 



(53) 



In this expression, the spectral dependance comes mainly from the effective polarizability 
and also from the far field reflection coefficients. One thus expects that the detected spectral 
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signal detected will follow polarizability variations that we have described in the preceding 
sections. 



In Fig. 12 we plot the signal detected in far field by a detector in the 45° direction from 
the vertical direction, in the case of a 100 nm-radius tungsten spherical particle heated at 
300 K and situated at different distances from the interface separating vacuum from SiC. 
When the particle-interface distance is 100 nm, the signal is peaked around SiC plasmon 
resonance, as is 100 nm-radius tungsten spherical particle polarizability. When the particle 
is retracted from the interface, the signal is reduced and broadens. Above a certain distance, 
the effective polarizability is nothing but the polarizability in vacuum. For tungsten, this 
polarizability has a rather flat dependance with the frequency. The signal spectral behaviour 
mainly comes from spectral variations of the reflection coefficient. One notes that the signal 
is more important in the frequency range where SiC is known to be highly reflective i.e. 
between 1.6x 10 14 rad s _1 and 1.8x 10 14 rad s _1 . A similar behaviour is observed for a 
micrometric tip except that the signal is broader at the minimum distance. This can easily 
be explained by inspecting effective polarizability of a 1 /zm-radius spherical tungsten tip 
situated at 1 /im of a SiC interface (tip in contact) where the noticed broadening is indeed 
observed. 

Note that for such situation, the detected signal can exhibit a peak very similar to the one 
observed in the energy density spectrum that is in the electromagnetic LDOS. However, this 
peak is the signature of the effective polarizability which exhibits a resonance at frequencies 
close to surface resonance. Even if such experiments do not probe EM LDOS, it can probe 
surface resonances if the tip is sufficiently close to the interface. It can then also be seen as 
a surface spectroscopy method. 



VII. RADIATIVE COOLING OF A DIPOLAR PARTICLE 

When a particle is heated, one can calculate the rate at which it will exchange energy 
with the outside. The spectral power lost by a particle at temperature T represented by its 
dipole when vacuum at null temperature is^ 

P(uj)du = -^Im[a(u)]e(u,T)<kj (54) 
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It has been shown in the past that when the particle is close to a surface it exchanges with 
it with a dependence of the polarizability and of the electromagnetic energy density at the 
particle position. This phenomenon is very similar to the Purcell effect when an atom or 
a molecule has its spontaneous emission rate that is modified when approached close to a 
surface. The cooling rate depends on the electromagnetic LDOS^Sl j n a similar way as the 
spontaneous emission rate depends on the EM LDOS. The cooling rate formula given in the 
paper of Mulet et alP^ should be corrected at close distances. Indeed, the interaction of the 
particle with the surface should be taken into account through the effective polarizability. 
The formula of the heat exchanged between the particle and the sample reads 

PMd W = ~9( W) T) E S{<4)St[G% B (r p ,rp)]du (55) 

i=x,y,z 



In Fig. 13 , we plot the heat exchanged between a SiC particle and a SiC substrate at 300 
K as a function of the tip-sample distance. We also show the spectrum of the exchanged 
flux. We compare, in each plot, the expression presented above and the one obtained by 
MuletpS. In the plot of the exchanged power as a function of the particle-surface distance, 
we observe that both expressions give similar results down to distances as low as 400 nm. 
At large distances, corrections to the polarizability are small . The exchanged heat comes 
mainly from the spectral region where a resonance occurs in the SiC particule, around 
935 cm -1 i.e. where e = —2. When one enters in the near field, two phenomena occur : 
electromagnetic field is more-strongly dominated by evanescent contributions, in particular 
those coming from the phonon polariton resonance occuring around the frequency where e=- 
1. Moreover, when the distance between the particle and the interface is reduced, corrections 
to the polarizability appear. In the present case, the particle resonance peak is weakened, 
which explains the fact that at 100 nm, the exact expression of the transfer gives a lower 
contribution than the one used by Mulet et alP^l. 



VIII. CONCLUSION 



We have shown how the interaction between a probing SNOM tip and the surface it 
scans modifies infrared electromagnetic near fields. We have shown that the particle dipole 
polarizabilities can be modified to account for this interaction, and we have added the effect 
of retardation to the one of the image dipole. We have analysed the relation between the 
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signal detected in the far-field by an apertureless SNOM and the excited fields, either due to 
external excitation of a plasmon or due to thermal emission. The technique performs a local 
spectroscopy of the surface and for some cases the signal spectrum can even be proportional 
to the EM LDOS, a fundamental quantity. We have also corrected the near-field radiative 
heat transfer formula for cases where there is a strong particle-sample interaction. In the 
future, this theory should be extended to tips with more elongated shapes^ such as conical 
ones, often used in experiment, in order to be applied to more realistic situations. 
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FIG. 11. Detected signal spectrum. The tip is constituted of a W spherical particle situated above 
Au. The signal is detected in a direction making an angle of 30° with the vertical (upper) Total 
signal and different contributions for a lOOnm radius tip at z =100nm above surface. (Lower) Total 
signal and different contributions for a 1 /xm radius tip at z =1 /um above surface. 
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FIG. 12. Signal at the detector which direction makes an angle = 45 PC with the vertical axis. 
The tip is a spherical tip of tungsten heated at T = 300 K and situated at various particle-SiC 
interface. Upper figure : signal in log-scale for Ru P = 100 nm and for three different tip-sample 
distance. Lower figure : signal in linear scale for Ru p = 1 /im and for separation distance equal to 
1 /j,m and 2.51 fun. Inset : signal for z = Ru P = 100 nm. 
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FIG. 13. (a) Heat transfer between a lOOnm radius SiC particle at 300 K above a SIC substrate 
at null temperature versus distance, (b) Power exchanged spectrum between a 100 nm radius SiC 
particle above a SiC substrate at zero temperature situated at 251 nm below the particle, (c) 
Power exchanged spectrum between a 100 nm radius SiC particle above a SiC substrate at zero 
temperature situated at 1 /im below the particle. 
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